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Abstract. A deceleration system consisting of staged parachute clusters and retro thrusters is optimized for the 
recovery of the first stage of a launch vehicle on sea. Optimal mass as well as reduction in speed by each parachute 
cluster and the retro thrusters is essential to minimize the inherent payload loss due to inclusion of additional 
systems. Three disciplines are involved in the study, namely parachute design, grain design and Three Degrees of 
Freedom (3-DOF) trajectory simulations. Parachute components are sized and their masses are estimated using a 
parachute design code. It computes the number of parachutes in the cluster, their sizes and opening loads for 
multiple reefing stages. Solid motor grain design is carried out, using high burn rate propellant, to provide high 
thrust to decelerate the launch vehicle stage to a near-zero descent rate at touchdown. A Multiobjective Mul¬ 
tidisciplinary Design Optimization (M^DO) problem has been formulated to minimize the mass of the deceler¬ 
ation system and minimize the touchdown speed of the recovered stage, subject to constraints on Maximum 
Expected Operating Pressure (MEOP), feasibility, etc. The optimization is carried out and the Pareto optimal 
front is obtained using an in-house multi-objective optimization algorithm. Attractor Anchored Multi-objective 
Evolutionary Algorithm (A^-MOEA). A total of twenty-hve design variables are considered including initial 
conditions for each deceleration stage, size of parachute cluster components for both drogue and main parachutes, 
and the size and shape of the solid motor grain for retro rockets. It is seen that the two objectives are conflicting. 
The Pareto optimal designs are discussed and the variation of design variables is presented. 

Keywords. Multidisciplinary design optimization (MDO); multiobjective multidisciplinary design 
optimization (M^DO); stage recovery; parachute cluster; solid motor grain design; three degrees of freedom (3- 
DOF) trajectory simulations; multi-objective evolutionary algorithm (MOEA); A^-MOEA. 


1. Introduction 

Recovery of the booster stage of launch vehicles has been 
studied extensively for reducing the cost of access to space. 
The high energy of the stage has to be dissipated before its 
touchdown on land or sea. This is carried out typically by 
either using the vehicle propulsion or employing aerody¬ 
namic decelerators. Eor final stage of the recovery, hover¬ 
ing, impact attenuators and even mid-air recovery 
techniques have been studied. A high Technology Readi¬ 
ness Level (TRL) stage recovery system named Parachute- 
Retro-Eloat (P-R-F) is proposed, which employs multi¬ 
stage parachute clusters followed by retro thrusters for 
landing with near-zero touchdown speed in sea on floats. 
The parachute system consists of the drogue and main 
parachute clusters. The drogue parachute decelerates the 
payload to a speed suitable for the deployment of the main 
parachute, which reduces the speed of the payload to the 
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desired terminal speed. The retro thrusters decrease the 
touchdown speed to the desired value before touching down 
in sea/land. It is essential to optimally size all the stages of 
the deceleration system together in order to reduce the total 
mass of the additional systems to ensure minimum payload 
loss. 

Yang et al{\] carried out optimization of reef parameters 
to minimize the maximum instantaneous opening load 
during the operation of a multi-reefed parachute using 
Genetic algorithm. Leonard et al [2] proposed a proba¬ 
bilistic algorithm for the optimization of drogue-to-main 
transition altitude for high altitude, low opening ballistic 
air-drop. Brooks [3] carried out grain configuration opti¬ 
mization based on ballistic design analysis for a star grain. 
Though parachute-propulsion soft-landing technique is well 
known [4], Multiobjective Multidisciplinary Design Opti¬ 
mization (M^DO) of the deceleration system has not been 
reported. 

This paper discusses the design of the deceleration sys¬ 
tem consisting of multi-stage parachute clusters (both 
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drogue parachute and main parachute) and the retro 
thrusters. The problem has been posed as a M^DO problem 
to minimize the mass of the deceleration system while 
simultaneously minimizing the impact or touchdown speed. 
The aim is to obtain the Pareto optimal set of deceleration 
system components along with their sizing, which opti¬ 
mizes these two objectives. 

In this paper we discuss the formulation of the M^DO 
problem and the salient features of the Pareto optimal 
designs. Section 2 presents the methodology. Section 3 
presents the optimal configurations and the variation of key 
design features. The conclusions are presented in Section 4. 


2. Methodology 

2.1 Optimization problem formulation 

The M^DO problem is posed in the Non-linear Program¬ 
ming Problem (NLP) form in Table 1, where Xdmgue, Xmam, 
Xretro refer to design variables of the drogue chute design, 
main chute design and the retro design disciplines, 
respectively. XcoupUng represents the set of coupling vari¬ 
ables. Xmin and Xmax refer to the lower and upper variable 
bounds, respectively. The first objective, /i, of the opti¬ 
mization problem is to minimize the sum of the masses of 
the drogue parachute, the main parachute and retro motors 
(Eq.(l)). The second objective, / 2 , is to minimize the 
touchdown speed (Eq. (2)). Eq. (3) and Eq. (4) quantify the 
feasibility or in-feasibility of a parachute system for a given 
set of design variables and fixed parameters. Occurrences 
like negative riser dimension, incompatible terminal state 
and reefing parameter, incompatible reefing fractions and 
disreef altitudes, and exceedance of maximum parachute 
diameter for a given maximum number of parachutes in a 
cluster constitute design in-feasibility. Eq. (5) checks for 
the existence of sustainable pressure inside the port volume 
for a given throat diameter. Eqs. (6) to (8) enforce 


Table 1. The M^DO problem 

minimize f\ . tn^rogue 3- niretro -f P (1 ) 

(^) fl ■■ ^touchdown 


X — {Xdrogue ^Xjrittin^Xretro ^Xcoupling} ^ \X,tiin t Xmcix\ 

subject to: 

• Xgrogue P h (3) 

g2 ■■ X^iin > 0 (4) 

83 • Xrctro ^ h (5) 

: 10’-M£OP> 0 (6) 

85 ■ De,max -De> 0 (7) 

86- rp- 20.0 > 0 (8) 


constraints on Maximum Expected Operating Pressure 
(MEOP), grain port radius and nozzle exit diameter. All the 
objectives and the constraints are normalized. 

The Individual Discipline Feasible (IDF) architecture has 
been employed in the present study to solve the Multidis¬ 
ciplinary Design Optimization (MDO) problem. A sche¬ 
matic of the M^DO architecture is shown in Fig. 1 . The set 
of twenty-five design variables (A) is optimized. X is 
divided into three subsets, Xdrogue, X,„am and Xretro, which 
are given to the three disciplinary modules: drogue para¬ 
chute, main parachute and grain design modules. Along 
with this a set of coupling variables, Xcoupling, is supplied by 
the optimizer to the disciplinary modules. The flow of 
design variables from the optimizer to the disciplinary 
analysis blocks is indicated by the blue arrows in Fig. 1. 
The constraint and objective function values are supplied 
back to the optimizer. This is indicated in green lines in the 
figure. The masses of the main parachute cluster and the 
retro motors are the coupling variables, which need to be 
fed back to the upstream modules. They are indicated 
below the modules along with dashed red lines. The esti¬ 
mated value of the coupling variables (Xcoupling = 
{mmam,est,mretro,e 5 :t}) IS Supplied by the Optimizer to the 
disciplinary analysis blocks, which is used for the designs 
of drogue and main parachutes. The difference between the 
estimated values and the computed values of Xcoupimg is 
used as penalty, given as follows: 


P = W2{m, 


manifest 


n ) -f W3{mr, 


^Xlobegnt ■ 


(9) 


The optimizer. Attractor Anchored Multi-objective Evolu¬ 
tionary Algorithm (A^-MOEA), employs real-coded 
Genetic Algorithm (GA) and cannot handle integer design 
variables at present; hence the number of lobes in a star- 
grain is represented as a continuous design variable in the 
optimizer, and a penalty equal to the square of difference 
between its value and its rounded off value {XNiobe,mtl is 
also added in Eq. (9). The penalty P is added to the first 



Figure 1. Architecture of the M’DO problem. 
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objective to ensure convergence of the coupling variables 
and the number of lobes. 

The first three modules in Fig. 1 contribute to the 
computation of the first objective function,/i, and all six 
constraints, g\-g6- The inputs from the grain design 
module go to the trajectory module, which evaluates the 
second objective function, / 2 . The terminal Equivalent 
Air Speed (EAS) of the drogue parachute, V/^, goes as 
an input to the main parachute design module as initial 
speed, V, 2 - Similarly, the terminal EAS of the main 
parachute, V/^, goes as an input to the trajectory simu¬ 
lation module as Vi^ = along with the retro thrust 

iT{t)) and propellant mass, mp{t). The trajectory code 
initiates the firing of the retro motor at an altitude of hi^, 
which is also a design variable. 

The design variable subsets Xdrogue and X^ain include 
suspension line length to nominal diameter ratio {Lo/Dq), 
forebody distance from parachute, reefing fraction, disreef 
altitude, terminal equivalent speed and terminal altitude. 
The deployment altitude for the drogue parachute is fixed; 
however, Xmain also includes parachute deployment altitude 
in addition to the design variables listed for the drogue 
chute. Xreiro iucludes retro ignition altitude and the grain 
design variables, namely outer diameter, length, port radius, 
slot tip height, slot base thickness, slot radius, ratio of port 
to nozzle throat areas, nozzle area ratio and number of lobes 
of grain. Initial state of drogue parachute deployment, body 
parameters, design factors and materials of parachute 
components, number of retro thrusters, nozzle half angle, 
bounds and other tolerances on various factors are fixed 
parameters for the problem. 


2.2 Parachute cluster design (Parachute Design 
and Analysis Code (ParDAC)) 

Parachute Design and Analysis Code (ParDAC) is an in- 
house code that designs parachute clusters to achieve a 
user-specified terminal condition from its state of 
deployment. When the parachute achieves terminal 
velocity, total weight of the parachute-body system is 
balanced by the total drag on the system. Hence, nominal 
diameter of the parachute is calculated via a force balance 
along the vertical direction. The drag coefficient of a 
parachute is dependent on a number of parameters like the 
shape of the canopy, rate of descent, the wake effect due 
to forebody, interference effect due to clustering of 
parachutes and so on. The experimental results presented 
by Knacke [5] are used to account parachute drag coef¬ 
ficient (Co„) as a function of rate of descent, the ratio of 
suspension line length to nominal diameter of parachute 
(Lo/Dq) and number of parachutes. Hence, drag coeffi¬ 
cient with respect to nominal area (5o) of the parachute is 
computed as 


Cdo ^^Orad ^ ^^^0 id- ^ N ( ^ O ) 

’ ^0 


where Cd„ ^ is Cdo at the user-defined rate of descent, 
J-Cd l, is a factor indicating the drag augmentation due to 

increase in ^ and J-Cdq a, is a factor quantifying drag 
reduction due to parachute clustering. 

Parachute opening force during each reefing stage is 
calculated using canopy loading (W/CdS) method pro¬ 
posed by Knacke [5], the drag area variation during para¬ 
chute inflation proposed by Yang et al{Y\ and the trajectory 
simulator described in Section 2.4. For a multi-reefed 
parachute, the drag area variation during its inflation to the 
first reefed configuration can be assumed to be linear with 
time, while the inflation following further disreefing is 
assumed to be exponential with time. Following [1], the 
variation of drag area (CoS) with time is taken as follows: 


CdS = < 


^^^t, {0<t<tp) 

tp 

(CdS)i, {tp<t<ty^) 


(CflS)i+[(CflS)2-(Cz,S)i] 

{tyy <t<tl) 


t-u 


y\ 


>/,i 


{CdS)„, (tn-l < t < f«-l+ty„) 

(CoS)„+[(CdS)„^,-(CoS),] , 

\ ^mn J 

, {tn—\ A ty^ A t A trt^ 

( 11 ) 

where tp, ty and t^ represent pre-inflation time, reefed time 
of a canopy reefing stage and Ailing time of a reefed stage, 
respectively. 

During inflation to the first reefed stage, opening force is 
given by Knacke [5] as follows: 


Fi 


(CdS), qCxX 
[(CdS)-(CdS),_,] q Cx Xi 


2, 3... 


( 12 ) 


where Cx and X refer to opening force coefficient at infinite 
mass (encountered in wind tunnel tests in which velocity 
reduction due to parachute is not achieved) and force 
reduction factor to incorporate the finite mass condition, 
respectively. In case of disreefing, opening force depends 
on fhe change in drag area during disreefing. The subscript i 
refers to the number of reefing stage for which opening load 
is calculated. CoS in Eq. (12) is obtained from reefing 
fractions that are design variables in optimization and the 
computed value of drag area of the designed parachute after 
complete inflation. Also, q at disreef altitudes is obtained 
from trajectory computations. 
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The maximum opening load (F,„ax) among the various F; 
is used for component sizing and mass estimation of vari¬ 
ous components. The design loads for canopy (fJ-icnpy), 
suspension lines {fisusp}’ radial lines {fJ-rad)^ ''^rit lines 
skirt band {fi.hand) and riser with Design 

Factor (DF) are as follows: 


Fmax ^ F)F (-fipy 


f^lcnpy 

^ max ^ cnpy 

n X {Dp) 

(13) 

f^susp 

^max ^ susp 

^susp 

(14) 

f^md 

F X DF 

^ max ^ susp 

^SLlSp 

(15) 

l^vline 

F X DF 

^ max susp 

^susp 

(16) 

l^sband 

O.OSFfyiQX ^ stand 

(17) 

f^riser 

— ^f^max ^ riser- 

(18) 


The number of material layers for the components is then 
decided using the material strength, followed by mass 
estimation. 


2.3 Ballistic Incremental Analysis Code (BINC) 

The Ballistic Incremental Analysis Code (BINC) evaluates 
the thrust time profile, pressure time profile and other 
characteristics of the solid motor grain. It is implemented as 
a C++ class to estimate the quantities required for the 
MDO study. For a given grain geometry, first the burning 
surface area of the grain and the port area are computed as a 
function of web thickness. The burn rate of the propellant is 
used to calculate the mass and energy flux into the control 
volume at a given time. Burn rate augmentation due to the 
flow properties is also accounted. Conservations of mass, 
momentum and energy along with experimentally obtained 
flame properties are employed to compute the one-dimen¬ 
sional flow inside solid motor. The pressure is iterated to 
convergence while accounting for the mass ejected through 
the nozzle. The module computes the vacuum thrust, pro¬ 
pellant mass and nozzle exit properties as a function of 
time. The overall flow of the BINC is shown in Fig. 2. 

The grain geometry evolution is estimated as a function 
of web thickness. This is accomplished using geometrical 
and analytical relations to get the exact perimeter and 
surface area of a particular grain. The grain geometry is 
decided based on the six design variables: grain outer 
dimeter, port_radius, slot tip height, slot base thickness, slot 
radius and number of lobes. This parameterization of the 
grain geometry results in hve different types of grains. 


namely slot-star with port, slot star, dimple star with port, 
dimpled star and cylindrical. In addition to the above six 
design variables, the grain length and the nozzle area ratio 
are the other propulsion related design variables. 

The computational domain, depicted in Fig. 3, consists of 
a number of sections divided along the length. One-di¬ 
mensional analysis is carried out and the gas state is esti¬ 
mated along the length of the grain. A database of the 
exposed surface area with constant increments in web 
thickness is calculated by BINC. It is used to estimate the 
surface area at any location and at any web thickness. The 
following assumptions are employed in order to estimate 
the pressure-thrust signature of the solid rocket motor 
(SRM): 

• Quasi steady 1-D analysis is carried out and no initial 
and final transients are computed. 

• Velocity at head end of the SRM is taken as zero. 

• No erosive burning is considered at the head end, as 
velocity is taken as zero. 

• The flow is assumed to be choked at the nozzle throat. 

• Ideal gas approximation for the combustion products is 
used. 

• Available experimental data has been employed for the 
throat degradation factor, burn rate augmentation 
factors and the flame properties. 

• Mukunda model [6] is used to account for erosive 
burning. 

• For the CG calculation, the section CG is assumed to 
be at the centre of the section. 

• Flame properties, the flame temperature, C*, Prandtl 
number, density, molecular weight, ratio of specific 
heats and absolute viscosity are only pressure 
dependent. 

• Thrust is assumed to be non-significant if the chamber 
pressure falls below 0.7 bar. 

During the analysis, the burn rate of the propellant at a 
given section is estimated using the Saint Roberts law as a 
function of pressure of the gases at that section. It is used to 
calculate the mass and energy flux fed into the control 
volume (section) at a given time. The burn rate is also used 
to calculate the incremental change in exposed surface area 
and port area at the next time step. The one-dimensional 
conservation of mass, momentum and energy equations 
along with experimentally obtained flame properties are 
employed to compute the special gradients of the pressure, 
temperature and density of gases. Burn rate augmentation 
due to the erosive burning is also accounted. The pressure is 
iterated to convergence to satisfy the 1-D conservation 
principals of the flow inside the SRM at a given time. 
Finally, the module computes the vacuum thrust, propellant 
mass and nozzle exit properties for a converged chamber 
pressure as a function of time. 
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Figure 2. Flow chart of the algorithm employed in BINC. 
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Figure 3. Schematic of a rocket motor with its grain divided into 
small elements. 


2.4 Three Degrees of Freedom (3-DOF) 

Trajectory Code 

The 3DOF equations are given in earth-centred rotating 
spherical co-ordinate system {r, 6, where r is the 
distance from the centre of earth, 9 is the longitude and f is 
the latitude (Fig. 4). The velocity is specified in terms of its 
magnitude (V), the flight path angle (y) and the heading 
angle (ip). The heading angle (i/t) is defined wifh respect to 
the local East in Velocity Aligned Reference Frame 
(VARF) (Fig. 5). First the three kinematic equations are 
listed, followed by the three dynamic equations: 

f = y sin(y) (19a) 



V cos(y) cos(i/^) — r cos(</)) (Oej 
r cos(0) 


(19b) 


V cos(y) sin(t/^) 
r 


(19c) 





Figure 4. Radius Aligned Reference Frame (RARF) 


V = -h ffl|.^rcos^(</))( sin(y) 

T 

— cos(y) tan(<^) sin(i/^)) H- 

' m 

L 

Vy = -cos(cr) — gcos(y) H-cos(y) 

m r 

+ IcoejV cos{(f>) cos(i/^) 

+ cos^((/)) (cos(y) + sin(y) tan((^) sin(i/^)) 

-C0S((T7') 

m 


(19d) 


(19e) 


L sm(cr) 

Vy/ = -^-cos(y) cos(i//) tan(0) 

mcos(y) r 

-F 2tt)£_/y(tan(y) cos((^) sin(i/^) — sin(</))) 

sin((/)) cos((/)) cos(i/^) (19f) 

cos(y) 

^ rL sin(crr) 
m cos(y) 

In these equations, D is the drag, L is the lift and T is the 
thrust. The thrust is resolved into two mutually orthogonal 
components - a component Ty along the positive velocity 
direction and the other component Ti perpendicular to it: 
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Ty = T cos 07 -, 

(20) 

Tp = T sin Op. 

(21) 


Here, dj is the angle between the thrust vector and the 
velocity vector. The lift bank angle is denoted by a and on 
similar lines, Gj denotes the thrust bank angle. The (Oej is 
the angular velocity of earth with respect to an Earth 
Centered Inertial (ECI) reference frame. The mass m is a 
variable that changes due to the change in propellant mass 
(m). The propellant mass flow rate can be computed using 
the specific impulse of the rocket, Isp, as follows: 

m = —. (22) 

Isp 

The trajectory code is implemented as a class in C++ 
language. The initial state vector, ?= {r,9,([>,V,y,\l/} is 
specified and these six equations are integrated using the 
RK-4 method. The spherical earth gravitational model has 
been used along with a standard atmosphere model with a 
non-rotating earth. 

2.5 Attractor Anchored Multi-objective 
Evolutionary Algorithm (A^-MOEA) 

In the current study we have used Attractor Anchored 
Multi-objective Evolutionary Algorithm (A^-MOEA) [7], a 
real-coded GA-based Multi-objective Evolutionary Algo¬ 
rithm (MOEA). Eor n/-dimensional objective space, the 
algorithm defines attractors in the space of first «/ — 1 
dimensions. In the selection operator of A^-MOEA, the 
solutions that are closest to these attractors are preferred 
over other solutions and are designated as Elite, which 
survive over GA generations. This ensures more uniform 
spread and better convergence of the Pareto front. 

3. Results 

The M^DO problem was solved for the recovery of spent 
rocket stage descending at an initial speed of 150 m/s at an 
altitude of 7 km. Three cases were considered, having spent 
stage masses of 4000, 8000 and 12000 kg. Population sizes 
of 140 and 250 were used for running A^-MOEA. Results 
of the optimization problems are presented in the form of 
Pareto fronts with Z-axis representing the first design 
objective, the T-axis representing the second design 
objective and the colour of the symbols representing the 
value of a third variable (Fig. 6). It can be seen that the two 
objectives are conflicting, i.e. reduction in terminal speed 
requires higher mass and vice versa. 

From the results (e.g., Fig. 6a) it can be observed that the 
Pareto front can be divided into three regions: Region Rl, 
where the terminal speed is high and the deceleration 


system employs only the drogue parachute; Region R2, 
where the terminal speeds are intermediate and the decel¬ 
eration is provided by both the drogue parachute and the 
main parachute clusters; Region R3, where the terminal 
speeds are low and retro propulsion is employed to reduce 
the terminal speed in addition to use of the drogue and main 
parachutes (Fig. 6b). 

This is because it is efficient to reduce the speed using 
parachutes for higher terminal speeds, where dynamic 
pressure is significant. As the mass of the parachute is 
inversely proportional to the square of the terminal speed, 
the terminal speed reduction using parachute becomes 
inefficient at low terminal speeds. This is seen from the 
results, where the terminal speed becomes relatively flat in 
the end of R2 region as a consequence of the large rise in 
main parachute diameter to reduce the terminal speed. 
However, the mass of propellant for a given AV is inde¬ 
pendent of the terminal speed. Therefore, at low terminal 
speeds, the optimizer configures a deceleration system that 
first employs parachutes followed by a retro motor to 
decelerate the stage. The slope of the Pareto front in region 
R3 becomes steeper as spent stage mass reduces, indicating 
diminishing advantage of retro system in impact speed 
reduction with increase in spent stage mass. The clear 
demarcation of regions Rl, R2 and R3 indicates that the 
optimization problem considered has variable dimension in 
design variable space. 

With increase in the spent stage mass, for a given ter¬ 
minal speed, the optimal deceleration system mass 
increases in all three regions, e.g., in region R3, larger 
retro motor is required to retard a heavier stage (Fig. 6b). 
Similarly, for a given mass of the deceleration system, the 
terminal speed increases with increase in spent stage mass. 
For a given spent stage mass, the size of drogue parachute 
increases with reduction in terminal speed in the region 
Rl (Fig. 6c). 

However, the drogue size is not sensitive to change in the 
terminal speed in the regions R2 and R3. Similarly, the size 
of main parachute rises with reduction in touchdown speed 
in region R2 (Fig. 6d). In region R3, the propulsion starts 
contributing while the main parachute reduces to a nearly 
constant size. The retro mass increases with the decrease in 
terminal speed. Transition between second and the third 
regions is influenced by the maximum allowable diameter 
of a parachute (Z)o,max =38 m) in cluster, which is a fixed 
parameter. The optimal design shifts from having no main 
parachute in Rl to systems with one, two or three main 
parachutes in R2 (Fig. 6a), and further to a single main 
parachute in R3. 

The decrease in terminal speed in regions Rl and R2 
leads to increase in the maximum opening loads on both 
drogue and main parachutes (Fig. 6e, f) due to increase in 
parachute size. The opening loads also increase signifi¬ 
cantly with increase in stage mass. Both the drogue and the 
main parachutes are reefed. 
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Figure 6. Results of deceleration system M^DO for recovery of 4000, 8000 and 12000 kg stages. 
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The reefing ratios and disreef altitudes for the drogue and 
the main parachutes have been optimized along with the 
deployment altitude of the main parachute and the ignition 
altitude for the retro motors. The optimal deployment 
altitude for the main parachute is shown in Fig. 6g and the 
optimal retro motor ignition altitude is shown in Fig. 6h. It 
is seen that when the terminal speeds are lower than wlO 
m/s (region R3) the deployment altitude of the main para¬ 
chute jumps to signihcantly higher altitudes along with a 
lower size of main parachute, as discussed earlier (Fig. 6d). 
It may be noted that lower size main parachutes, in the 
region R2, are also deployed at higher altitudes. The opti¬ 
mal ignition altitude for the retro motor is very small in all 
cases, 8-18 m in the present case, and does not vary sig¬ 
nificantly with decrease in terminal speed. The optimal 
ignition altitude increases with increase in the spent stage 
mass, as expected. 


4. Conclusions 

The design of a deceleration system for stage recovery 
has been posed as an M^DO problem. The deceleration 
system two-stage parachutes (with reefing) and a retro 
motor. Two parachute design modules (for design of 
drogue and main parachute clusters), a retro motor design 
module and a trajectory module have been implemented 
for this multidisciplinary design optimization problem. 
Two conflicting objectives, namely the mass of the 
deceleration system and the touchdown speed, are mini¬ 
mized together, subject to constraints in parachute design 
and retro motor design while considering twenty-five 
design variables. 

The optimal design set divides into three distinct subsets. 
First subset: the optimal deceleration system for touchdown 
speeds larger than k, 25 m/s consists only of drogue para¬ 
chutes; second subset: for touchdown speeds between ?»10 
and 25 m/s, the optimal design solution consists of only 
drogue and main parachutes and third subset: for terminal 
speeds below «10 m/s, it is optimal to have drogue para¬ 
chutes with retro motors and less number of main para¬ 
chutes. Another interesting observation is that the size of 
drogue parachutes is insensitive to the desired touchdown- 
speed for the second and third subsets. Therefore, the pre¬ 
sent M^DO study is able to discover broad design rules for 
this multidisciplinary design problem. 

Thus depending on the impact speed of interest, the 
recovery system design optimization can be simplified by 
choosing the optimal deceleration system constituents 
which would reduce the number of design variables for 
regions R1 and R2. 

The masses of drogue parachute, main parachute and 
retro motor systems increase with decrease in terminal 
speed in the regions where these systems dominate the 
deceleration of the stage. The optimal variation of various 


paramaeters related to the parachute and retro motor 
deceleration phases have been discussed. 


Nomenclature 


Acronyms 

3-DOF 

Three Degrees of Freedom 

A^-MOEA 

Attractor Anchored Multi-objective 

BINC 

Evolutionary Algorithm 

Ballistic Incremental Analysis Code 

CG 

Center of Gravity 

DE 

Design Factor 

EAS 

Equivalent Air Speed 

ECI 

Earth Centered Inertial 

GA 

Genetic Algorithm 

IDE 

Individual Discipline Eeasible 

M^DO 

Multiobjective Multidisciplinary Design 

MDO 

Optimization 

Multidisciplinary Design Optimization 

MEOP 

Maximum Expected Operating Pressure 

MOEA 

Multi-objective Evolutionary Algorithm 

NLP 

Non-linear Programming Problem 

P-R-F 

Parachute-Retro-Float 

ParDAC 

Parachute Design and Analysis Code 

TRL 

Technology Readiness Level 

VARF 

Velocity Aligned Reference Frame 

Greek symbols 


(/) Elevation Angle, Latitude (rad) 
y Flight Path Angle (rad) 

6 Longitude (rad) 

0^ Angle between thrust vector and the velocity vector 
(deg) 

(o Angular Velocity (rad/s) 

xj/ Heading Angle (rad) 

O' Bank Angle for lift (deg) 

cr^ Bank Angle for thrust (deg) 

t] Opening profile shape exponent 

1.1 Design load 

fj Design load per unit length 

Other symbols 

CdS Drag area 

Cdo Drag coefficient defined with respect to area 

corresponding to nominal diameter 
Coo KOD ^Do ^ particular rate of descent 
Cx Opening force coefficient at inhnite mass 

(dimensionless) 

D Diameter 

Dq Nominal diameter 

Do^nuix Maximum allowable diameter of a parachute in 
a cluster 

De Nozzle exit diameter 
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lobe ,int 


TC, 


n 


V 

Ar, 


^drogue 


De 

Dp 

F 

h 


hp 

L 

Le 


m 

rh 

C* 


blsusp 

p 

q 

rp 

r 

So 

s 

T 

t 


ly 

Tl 

Tv 

V 

W 


w 

X 

X 


Penalty for non-integer value of the number of 
lobes 

Factor indicating drag augmentation due to 
increase in It is 1 when ^ — 1 
Factor indicating drag reduction due to 
parachute clustering. It is one when number of 
parachutes in cluster is 1 
Equivalent Air Speed 

Indicates feasiblity of sustainable combustion 
inside the retro motors 
Indicates feasibility of a parachute design 
Indicates feasibility of a parachute design 
Maximum nozzle exit diameter 
Projected diameter of parachute 
Opening load 
Height, m 

Specihc Impulse, Ns/kg 
Length 

Length of suspension line 
Mass, kg 

Mass flow rate of propellant (kg/s) 
Characteristic velocity of the propellant 
Number of objective functions 
Number of suspension lines 
Penalty 

Dynamic Pressure 
Grain port radius 
Radius (m) 

Reference area corresponding to Do- So = ^ 
The state vector s = {r, 6, (f), F, y, ij/} 

Thrust, N 
time, s 

Pilling time of a reefed stage 
Pre-inflation time 

Reefed time of a canopy reehng stage 
Component of thrust along vehicle lift (N) 
Component of thrust along vehicle velocity (N) 
Speed (m/s) 

Weight carried by canopy 
weight 

Design Variable 

Porce reduction factor (dimensionless) 


Subscripts 


1 

For drogue parachute 

2 

For main parachute 

3 

For retro bring 

cnpy 

Canopy 

coupling 

coupling 

drogue 

Drogue parachute 

est 

Estimated 

f 

Terminal 

i 

Initial 

main 

main parachute 

max 

Maximum value 

P 

Propellant 

rad 

Radial line 

retro 

retro rocket 

riser 

Riser 

stand 

Skirt band 

susp 

Suspension line 

vline 

Vent line 
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